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The cross section of the process e~^e~ — > 3(7r+7r~) has been measured 
using 22 pb~^ of integrated luminosity collected with the CMD-3 detector at 
the VEPP-2000 e~^e~ collider in the cm. energy range 1.5 - 2.0 GeV. The 
measured cross section exhibits a sharp drop near the pp threshold. A first 
study of dynamics of six-pion production has been performed. 



1. Introduction 



Production of six pions in e~^e~ annihilation was studied at DM2 
and with much larger effective integrated luminosity at BaBar |2|, using 
Initial-State Radiation (ISR) events. The DM2 experiment observed a "dip" 
in the cross section at about 1.9 GeV, confirmed later by the FOCUS Col- 
laboration in the photoproduction [s], 0] and by the BaBar Collaboration, 



where this structure was also observed in the 2('7r+'7r~)'7r'^7r° final state [2|. 
The origin of the "dip" remains unclear, but the most popular explanation 
suggests a presence of the under-threshold proton-antiproton (pp) resonance. 
This hypothesis is supported by the fast increase of the pp form factor to 
the threshold, recently confirmed by the high-statistics BaBar study js], and 
discussed in many theoretical papers (see, e.g., Ref. {gI). Even earlier, a nar- 
row structure near the proton-antiproton threshold has been also observed 
in the total cross section of e"''e~ annihilation into hadrons in the FENICE 
experiment 

The e~^e~ — ?■ 3(7r"^7r~) cross section is also used in the calculation of the 
hadronic contribution to the muon anomalous magnetic moment Q]. The 
detailed study of the production dynamics can further improve the accuracy 
of these calculations and help in explaining the cross section anomaly. 

In this paper we report the analysis of the data sample based on 33 pb~^ 
of integrated luminosity collected at the CMD-3 detector in the 1.0-2.0 GeV 
center-of-mass energy range. We observe only a few candidate events below 
1.5 GeV. Since their number is consistent with background, we present our 
results for the 1.5-2.0 GeV center-of-mass energy range, corresponding to 
22 pb~^ of integrated luminosity. These data were collected in three energy 
scans performed at the VEPP-2000 collider jol. 
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The general purpose detector CMD-3 has been described in detail else 



where 10|. Its tracking system consists of a cylindrical drift chamber (DC) 11| 
and double-layer multiwire proportional Z-chamber, both also used for a trig- 
ger, and both inside a thin (0.2 Xq) superconducting solenoid with a field of 
1.3 T. The liquid xenon (LXe) barrel calorimeter with 5.4 Xn thickness has 
fine electrode structure, providing good spatial resolution [l2|, and shares 
the cryostat vacuum volume with the superconducting solenoid. The barrel 
Csl crystal calorimeter with thickness of 8.1 Xq is placed outside the LXe 
calorimeter, and the end-cap BGO calorimeter with a thickness of 13.4 Xq is 
placed inside the solenoid [l3|. The luminosity is measured using events of 
Bhabha scattering at large angles [l4|. 



2. Selection of e+e — >• 3(7r+7r ) events 



Candidates for the process under study are required to have five and more 
charged-particle tracks with the following "good" track definition: 

• A track contains more than five hits in the DC. 

• A track momentum is larger than 40 MeV/c. 

• A minimum distance from a track to the beam axis in the transverse 
plane is less than 0.5 cm. 

• A minimum distance from a track to the center of the interaction region 
along the beam axis Z is less than 10 cm. 

• A track has a polar angle large enough to cross half of the DC radius. 

The number of events with seven or more selected tracks is found to be 
less than 1%. Reconstructed momenta and angles of the tracks for six-track 
and five-track events were used for further selection. 

For six- or five-track candidates we calculate the total energy and total 
momentum assuming all tracks to be pions: 

5,6 5,6 

Etot = J2 VpF+^ , Ptot = I ^pi|. 

i=l i=l 
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Figure 1: (a) Scatter plot of the difference (AE) between the total energy and cm. 
energy versus total momentum for six-track events. The line shows the boundary of the 
applied selection; (b) Projection plot of (a) after selection. The histogram shows the 
normalised MC-simulated distribution; (c) Scatter plot of difference {AE) of total energy 
and cm. energy versus total momentum for five-track events; (d) Difference between 
the total energy of five-tracks plus missing track energy and cm. energy (points). The 
histogram shows the distribution for the MC simulated background events (see text). 



Figure [U^a) shows a scatter plot of the difference between the total energy 
and cm. energy AE=Etot-Ec.m. versus total momentum for six-track can- 
didates. The histograms combine events from three highest energy points. A 
clear signal of six-pion events is seen as a cluster of dots near zero. Events 
with a radiative photon have non-zero total momentum and total energy 
which is always smaller than the nominal one. A momentum of any pion 
incorrectly reconstructed due to interaction with detector material or DC 
resolution leads to momentum-energy correlated "tails" in both directions. 
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We select events with total momentum less than 150 MeV/c and show the 
difference AE in Fig. [^b). The experimental points are in good agreement 
with the corresponding Monte Carlo (MC) simulated distribution shown by 
the histogram. We require -200< AE<100 MeV to determine the number of 
six-pion events. Six-track events have practically no background: we estimate 
it from MC simulation of the major background processes 2(7r"'"7r~7r°) and 
2(7r+7r^)7r° (one of the photon from the 7r° decay converts to a e~^e~ pair at 
the vacuum pipe), and found a contribution of less than 1%. We use this 
value as an estimate of the corresponding systematic uncertainty. 

To determine the number of six-pion events with one missing track, a 
sample with five selected tracks is used. A track can be lost if it flies at 
small polar angles outside the efficient DC region, decays in flight, due to 
incorrect reconstruction, nuclear interactions or by overlapping with another 
track. Figure Wic) shows a scatter plot of the difference AE between the 
total energy and cm. energy versus total momentum for five-track events. 
Six-pion candidates in the five-track sample have energy deficit correlated 
with the total momentum. This sample has some admixture of background 
events from multihadron processes mentioned above with photons from the 
7r° decays. We apply an additional requirement on the "neutral" (not asso- 
ciated with charged tracks) energy Eneutrai in the calorimeter to be less than 
300 MeV. This requirement reduces the background by a factor of two and 
removes less than 2% of signal events estimated using MC simulation. 

The direction and momentum of a missing pion can be calculated assum- 
ing a six-pion final state. We add energy of a missing pion to the energy 
of five detected pions and show the difference AE in Figure [U^d) by points. 
A corresponding background distribution from the MC simulation of the 
2(7r+7r^7r°) and 2(71^71^)71:^ events is shown in Figure [U^d) by the histogram: 
background events contribute less or about 10% to the signal region after 
applying a requirement Eneutrai <300 MeV. 

To obtain the number of six-pion events from the five-track sample, we fit 
the distribution shown in Fig[T]^d) with a sum of functions describing a signal 
peak and background. The signal line shape is taken from the MC simulation 
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of the six-pion process and is well described by a sum of two Gaussian distri- 
butions. The photon emission by initial electrons and positrons is taken into 
account in the MC simulation and gives a small asymmetry observed in the 
distributions of Figs. [1] (b,d). We describe this asymmetry by an admixture 
of a third Gaussian function. All parameter ratios of the signal function are 
fixed except for the number of events and main Gaussian resolution. The 
third-order polynomial is used to describe the background distribution. 

To estimate a systematic uncertainty of the background subtraction pro- 
cedure, we compare the MC simulated background distribution with the 
experimental events with an Eneutrai > 300 MeV requirement, and found 
reasonable agreement with the histogram shown in Fig[T]^d). A variation of 
the polynomial fit parameters for the experimental and MC simulated back- 
ground distributions leads to about 3% uncertainty on the number of signal 
events. 

We found 2887 six-track events and 5069 five-track events corresponding 
to the process e+e^ — )■ 3(7r+7r^). The numbers of six- {Nq^^) and five-track 
{N^T^) events determined at each energy point are listed in Table [U 



3. First study of the production dynamics 

To obtain a detection efficiency, we simulate six-pion production in a 
primary generator, pass simulated events through the CMD-3 detector using 



the GEANT4 15| package, and reconstruct them with the same reconstruc- 
tion software as experimental data. In our experiment, the acceptance of 
the DC for the charged tracks is not 100%, and the detection efficiency de- 
pends on the production dynamics of six pions. The dynamics of the process 
e~^e~ — )■ 3(7r+7r^) was not studied previously in detail. The BaBar Collabo- 
ration ^ reported the observation of only one p(770) from all vr+vr" invariant 
mass combinations and no structures in any other (three-, four-pion) invari- 
ant mass combinations. 

We investigate a few production mechanisms, and compare simulated 
angular and invariant mass distributions with those in data. All studied 
distributions strongly contradict to a phase space model, which assumes all 
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pions to be completely independent. We exclude the phase space model 
from further consideration. In this paper we illustrate our study with three 
models, all with one p(770) per event. To conserve the initial state quantum 
numbers, six pions must have J^^ = 1 . 




Figure 2: Cosines of the relative angle of two pions with opposite-sign charge (a), and 
of two pions with same-sign charge (b) for experimental events (dots) and MC simulation 
for p(1420)7r+7r- ^ ai(1260)7r7r+7r- (dotted histogram), p(770)/o(1370) (solid histogram) 
and p(770)/2(1270) (dashed histogram). 



In the model #1 we use the following decay chain: e'^e — )■ p(1420)(7r+7r )s-wave 
ai(1260)±7r^7r+7r- p(770)°2(7r+7r-) 3(7r+7r-). This model uses dom- 



inant decays p(1420)° ai(1260)±7rT and ai(1260)± p(770)°7r± [16|, 
and naturally includes the ai(1260)^ — p(770)^7r° decay to describe the 
e"'"e~ — >■ 2(7r+7r^7r°) process with one charged p(770) [3]. We use PDG val- 



ues 



17| for the resonance parameters and the model allows to introduce a 



form factor in each decay vertex. 

Another studied model (#2) was simpler: it includes the production of 
one p(770)° and four pions in S-wave. We try two options: the four pions 
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are distributed according to the phase space or forming a scalar resonances 
/o(1370) or /o(1500). 

And finally, the model (#3) assumes e~^e~ — )■ p(770)/2(1270) with a ten- 
sor /2 resonance in the four-pion final state. 

MC simulation should reproduce experimental angular distributions of 
the pions to obtain correct detection efficiency. Figure [2] shows (by points) 
the cosines of open angles between pions for opposite-sign (a) and same-sign 
(b) pion pairs for data. 

We compare distributions of Fig. [2] with the MC simulated distributions 
for the model #1 (dotted histogram), model #2 (solid histogram) and model 
#3 (dashed histrogram), and the best agreement was found with the model 
#2. 

Note, that variation of the resonance parameters in the models does not 
significantly affect these angular distributions. For example, model 7^2 with 
production of one p(770)° exhibits the same angular distributions both in 
the case, when the remaining four pions are distributed according to phase 
space or form a scalar resonance (/o(1370) or /o(1500)). 



Figure [3](a) presents the polar angle (6',^) distribution for six-pion events 
with all detected tracks. The requirement for a track to cross half of the 
DC radius effectively determines a cut on this parameter. The result of 
the MC simulation in model #2, presented by the histogram, well describes 
the observed distribution. Figure El^b) presents the polar angle distribution 
for five detected tracks (circles for data, the solid histogram for the MC 
simulation) after background subtraction. The polar angle distribution for 
the missing track is shown by squares (data) and the dashed histogram (MC). 
With our "effective" DC acceptance we have almost two times more six-pion 
events with one missing track than events with all tracks detected. 

We calculate invariant masses for the combinations of two, four (total 
charge zero), and three (total charge ±1) pions for the different cm. energies 
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Figure 3: (a) Polar angle distribution for six-pion events with six detected tracks for 
data (points) and MC simulation (histogram); (b) Polar angle distribution for six-pion 
events with five detected tracks for data (circles) and MC simulation (solid histogram). 
The polar angle distribution for a missing track is shown by squares (data) and the dashed 
histogram (MC simulation). 

and show them in Fig. HI We compare the obtained distributions with model 
7^2 (p47r), and observe good agreement with experiment at cm. energies 1600 
MeV and 2000 MeV, if four pions are distributed according to phase space 
(sohd histogram). But at the cm. energy of 1800 MeV the experimental data 
are better decsribed by the same model with four pions forming /o(1370). 
Note that invariant mass distributions for models #1 and #3 do not describe 
data in any mass interval, but some admixture of these channels cannot be 
excluded. 

From the study of the mass distributions in Fig. H] we conclude that pro- 
duction djTiamics of six charged pions changes in the relatively narrow energy 
region (1700-1900 MeV). This phenomenon demands a further investigation. 



Ec.m. = 1.6 GeV 




Figure 4: Experimental invariant mass distributions (from left to right) for two, three and 
four pions for (top to bottom) 1600, 1800 and 2000 MeV cm. energies in comparison 
with simulation of one p{770)'^ with the remaining four pions in S-wave, and distributed 
according to the phase space (solid histogram) or form a scalar resonance /o(1370) (dashed 
histogram). 
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4. Detection efficiency 

We calculate the detection efficiency from the MC simulated events 
as a ratio of events after selections described in Sec. |2] to the total number 
of generated events. With the limited DC acceptance, incorrect simulation 
of the pion angular distribution leads to a systematic error in the efficiency 
calculation and thus in the cross section measurement. 




Figure 5: (a) Ratio of events with five and six detected tracks for data (points with, 
errors) and MC simulation for model #1 (dotted line), model #2 (solid line) and model 
#3 (dashed hue); (b) Double ratio Rif''/Rll^*^ versus energy. The line shows a fit with 
a constant. 



In the five-track sample, about 15-17% of events have a missing track due 
to the DC reconstruction inefficiency, well reproduced by the MC simulation. 
The remaining events migrate from the six- to the five-track sample due to 
the limited DC acceptance (see Fig. |3]). It makes the ratio R^q = N^t^/Nq^^ 
very sensitive to the pion angular distribution, and we study it to validate 
the model used for the efficiency calculation. 

Figure |5] (a) shows the i?56 ratio versus energy for data (points with errors) 
and for three models, discussed in Sec. El The experimental average value 
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TDdata 
^56 



MC#2 
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56 



1.76 for the model 
= 1.92, dotted line) 



1.74±0.03 is in good agreement with i?5g 
#2 (solid line), but inconsistent with model #1 {R 
and model #3 (i?^'"*^ = 1.30, dashed line). A "naive" phase space model 
for the six-pion production (all tracks uncorrelated) gives -R^*" =2.1. 

To estimate a model-dependent systematic error, we compare the exper- 
imental number of six- and five-track events after normalisation to the MC 
simulated acceptance. We calculate a double ratio Rff"'/R'^l^*'^ for each 
energy point for the model #2, and show it in Fig. O (b). The average value 
0.984±0.018 {x^ / n.d.f =bQ is in good agreement with the prediction of 
model #2 in the studied energy interval, so that a maximum systematic 
deviation from unity does not exceed 3.4%. However, a relatively large 
value can be an indication of the additional systematic uncertainty, and we 
conservatively take 4% as an estimate of a systematic error on the detection 



efficiency using Jx^ /n.d.f as a scale factor. 
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Figure 6: (a) Efficiency calculated from the MC simulation for six-track events (squares) 
and for a sum of five- and six-track events (circles); (b) Radiative correction. 



The detection efficiency thus obtained with model #2 is shown in Fig. |6t^a) 
for events with six detected tracks (squares) and for a sum of five- and six- 
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track events (circles), icreasing efficiency by factor 2.5. Note that if a sum of 
six- and five-track events (Nqt^ + N^t^) is taken for the detection efficiency cal- 
culation, the data-MC inconsistencies in the description of the DC inefficiency 
and (partly) in the model-dependent angular distributions are significantly 
reduced. 

5. Cross Section Calculation 

At each energy the cross section is calculated as 

L-e-(l + 5)' 

where L is the integrated luminosity for this energy point, e is the detec- 
tion efficiency (Fig. EJ^a)), and (1 + 5) is the radiative correction calculated 



according to [18| and shown in Fig. [6] (b). The energy dependence of the 
radiative correction reflects a sharp dip in the cross section. To calculate the 
correction we use BaBar data [2] as a ffist approximation and then use our 
cross section data for iterations. 

The integrated luminosity, the number of six and five-track events, detec- 
tion efficiency, radiative correction and obtained cross section for each energy 
point are listed in Table [1] 



6. Systematic errors 

The following sources of systematic uncertainties are considered. 

• The model dependence of the acceptance is determined using the an- 
gular distributions, which are specific for each particular model. As 
shown in Sec. HI a model with one p(770) and remaining pions in S- 
wave (phase space or /o(1370)) gives good overall agreement with the 
observed angular distributions. Using the ratio of six- and five-track 
events we estimate a systematic uncertainty on the detection efficiency 
as 4%. 
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Figure 7: The e+e" — > 3(7r+7r~) cross section measured with the CMD-3 detector at 
VEPP-2000 (dots). The results of the BaBar measurement Q are shown by open circles. 
The line shows the pp threshold. 



Since only one charged track is sufficient for a trigger (99-98% effi- 
ciency), we assume that for the multi-track events, considered in this 
analysis, the trigger inefficiency gives a negligible contribution to the 
systematic error. 

A systematic error due to the selection criteria is studied by varying 
the cuts described previously and doesn't exceed 3%. 

The uncertainty on the determination of the integrated luminosity 
comes from the selection criteria of Bhabha events, radiative correc- 



tions and calibrations of DC and Csl and does not exceed 2% 14 1 



The admixture of the background events not subtracted from the six- 
track sample is estimated as 1%. 
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• The accuracy of background subtraction for five-track events is stud- 
ied by the variation of functions used for a background description in 
Fig. dl^d) and is estimated as 3%. 

• A possible uncertainty on the beam energy is studied using the mo- 
mentum distribution of Bhabha events and total energy of four-pion 
events. The uncertainty at the level of 5 ■ 10~^ is not excluded and 
because of the cross section variation it can result in a 1% change of 
the cross section. 

• A radiative correction uncertainty is estimated as about 1% mainly 
due to the uncertainty on the maximum allowed energy of the emitted 
photon, as well as from the uncertainty on the cross section. 

The above systematic uncertainties summed in quadrature give an overall 
systematic error of about 6%. 

The obtained cross section is in overall agreement with the results of the 
most precise measurement performed by the BaBar Collaboration ^ shown 
in Fig. [7] by open circles. 

Conclusion 

The total cross section of the process e~^e~ — t- 3(7r+7r^) has been mea- 
sured using 22 pb~^ of integrated luminosity collected by the CMD-3 detector 
at the VEPP-2000 e^e~ collider in the 1.5-2.0 GeV cm. energy range. The 
five- and six-track events are used to estimate the model-dependent uncer- 
tainty in the acceptance calculation. From our study we can conclude that 
the observed production mechanism can be described by the production of 
one p(770) with four remaining pions in S-wave and distributed according to 
phase space. We also observe that the production dynamics changes in the 
1700-1900 MeV c.m.energy range and demands further investigation. A de- 
tailed analysis of the production dynamics will be performed in the combined 
analysis of the processes e+e~ — )■ 3(7r+7r~) and e~^e~ — )■ 2{n~^n~)27i^. 
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The measured cross section is in good agreement with all previous ex- 
periments in the energy range studied, and exhibits a sharp dip near the pp 
threshold. 
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Table 1: Luminosity, Number of events, Detection efficiency, Rad. correction and Cross 
section for each cm. energy point. Horizontal lines separate three energy scans. 
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